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Abstract — Nanofiltration is currently applied in many 
industrial processes. The separation efficiency of nanofiltration 
systems is related to complex phenomena occurring at 
membrane surface and within the nanopores. The nature of 
these phenomena is still a subject of debate and there is a real 
need to better reproduce these phenomena through simple and 
accurate predictive models. In this paper, interfacial and 
dielectric properties of two commercial nanofiltration 
membranes are investigated with the modeling of the 
permeation of ions typically found in seawater. The membrane 
charge density was estimated using zeta potential measurements 
and the dielectric exclusion was represented by the Born model. 
The predictions of rejection and permeate flux were in good 
agreement with experimental results when the dielectric effect 
was considered, indicating that the calculation of membrane 
charge with zeta potential data is appropriate. Based on 
simulation results, dielectric constants inside nanopores were 
calculated and results show that the ion solvation model is 
appropriate for these membranes. 

Index Terms — Dielectric exclusion, Ionic rejection, 
Nanofiltration, Zeta potential. 


I. Introduction 

Nanofiltration (NF) is a pressure-driven membrane 
separation process with characteristics between those of 
reverse osmosis and ultrafiltration and is currently applied in 
many industrial processes such as desalination [1]. The 
separation efficiency of nanofiltration systems is related to a 
complex mechanism including steric, dielectric and 
electrostatic partitioning effects between membrane and 
solutions [2], [3]. 

During the last two decades, the prediction of membrane 
performance has been a relevant area of research [1]. There is 
an increasing need for developing of model-based tools to 
design new membrane systems or to optimize existing 
membrane installations. These models should predict fluxes 
and rejections as a function of transmembrane pressure for a 
given membrane system. Also, they should be able to 
determine the membrane properties necessary to attain a 
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desired retention and be a realistic predictive tool with a 
limited number of experiments [3], [4]. 

The most widely and successfully adopted NF predictive 
models are based on the extended Nernst - Planck (ENP) 
equation to describe the mass transfer across the membrane 
[3], [5], [6]. This model considers the three important 
mechanisms of ionic transport in membranes: (a) diffusion, 
(b) electromigration as a result of concentration and electrical 
potential gradients and (c) convection caused by the pressure 
difference across the membrane [7]. One of the most studied 
models is the Donnan-Steric Pore Model (DSPM) [8], [9]. 
This model describes the transport of ions in terms of an 
effective membrane thickness (Ax), a membrane charge 
density (X d ) and an effective pore radius (r p ) [5], [6], [8], [9]. 
It also takes into account the effects of hindrance to diffusion 
and convection within the pore and the equilibrium 
partitioning due to a combination of Donnan and sieving 
mechanisms at membrane / solution interfaces [6]. Although 
this model has been reported to successfully describe simple 
systems such as those constituted by organic molecules, it has 
not been very successful for multivalent cations. To improve 
the prediction capability, some modifications for DSPM 
model have been suggested by Bowen and Welfoot [9] such 
as the incorporation of dielectric constant variations between 
bulk and pore solutions, which has shown better prediction of 
divalent ions rejection. Bandini and Vezzani [10] proposed a 
more general model, called Donnan-Steric Pore Model & 
Dielectric Exclusion (DSPM&DE), which is basically an 
extension of the DSPM model, in which the primary effect of 
the DE is considered as the most relevant in determining ion 
partitioning, together with steric hindrance and Donnan 
equilibrium. 

The membrane charge density is obtained by fitting 
rejection data in the DSPM and DSPM&DE models, being an 
empirical function related to the feed electrolyte 
concentration in terms of a Freundlich isotherm [10], [11]. 
This model is independent of the electrolyte type and does not 
consider any pH effect. It has been demonstrated to be 
appropriate in the case of single salts and multicomponent 
mixtures for some membranes, but it failed in some other 
cases [2]. Hence, it has been suggested that the membrane 
charge density is related to zeta-potential data by measuring 
the streaming potential of nanofiltration membranes, 
considering the influence of the ionic strength and pH. Other 
possibility is to develop physico-chemical models to describe 
the mechanism of charge formation, considering dissociations 
of functional groups [2], [12]. 

In this study, a model based on the DSPM models 
equations is used to predict the rejection of various ions 
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typically present in the seawater permeating through two 
different nanofiltration membranes. Zeta potential is 
incorporated in the model to estimate the membrane effective 
charge density for each salt solution studied. In the model the 
osmotic effect is considered by using the van’t Hoff equation, 
while the dielectric exclusion is expressed in terms of the 
solvation energy barrier. The Born term and a variation of 
viscosity within the pore are also taken into account. The 
dielectric constant inside the pore (e p ) is adjusted 
experimentally and the applicability of the single layer 
oriented molecules as the dominant dielectric exclusion 
mechanism is evaluated for both membranes. Finally, the 
simulation results using the DSPM-based models considering 
and not considering the dielectric exclusion mechanism are 
discussed and compared to the experimental data in order to 
validate the modified models for single salts solutions. 

II. THEORY 

A. Transport Equations 

The Donnan-Steric Pore Model (DSPM), originally 
developed by Bowen et al. [8], describes the ionic transport 
across a membrane by the extended Nerst-Planck equation 
(ENP) [9]. Steric effects are caused by the difference between 
membrane pore radius and the ion Stokes radius, while 
electrical (Donnan) effect is the result of charge distribution 
in the membrane and in bulk solution. The combined effect 
determines the selective ions transport trough the membrane. 

The ionic transport involves convection, diffusion and 
electric forces resulting in an ionic flux j t through the 
membrane that is represented by ENP [9]: 


Ji =K^ Ci U—^ 


c i D w „ K ,,i d/u, 


RT 


dx 


( 1 ) 


where K i c and K l d are the hindrance factors for convection 
and diffusion of ion i, respectively, D i oo is the bulk diffusion 
coefficient of ion i (m 2 s -1 ), p t is the electrochemical potential 
of ion i (J mol' 1 ), R is the universal gas constant (J mol' 1 K' 1 ) 
and T is the temperature (K). 


K i c = (1 + 2Z. - Z. 2 ) 

X (1 + 0.054/1,. - 0.988/1, 2 + 0.44 U,. 3 ) 

K id = 1 - 2.3/1, +1.154/1,. 2 +0.224/1, 3 

The mass transfer is steady state, so that the mass 
accumulation rate is null. Thus the molar fluxes through 
boundaries must be equal to fluxes inside membrane [15]. The 
ionic flux of component i correlates the molar flux of 
component i at permeate (Q p ) to solvent average velocity 
inside membrane [6], [9], [13], [16]: 


(4) 

(5) 


Substituting (2) and (6) into (1) yields the concentration 
gradient within the pore: 


dc t 

dx 
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where c t is the concentration within the pore (mol m' 3 ), D ip is 
the pore diffusion coefficient (mV 1 ), V, is the partial molar 
volume of i, Zi is the valence of ion i, F is the Faraday constant 
(96487 C mol' 1 ), y/ is the electrical potential within the pore 
(V) and C it f and C i>p are the measured solute concentration in 
the feed and permeate sides, respectively. 

The electroneutrality conditions at feed side, inside 
membrane pore and at permeate side are respectively: 


2>A/ =° 



2>A,=o 


i =1 


(8) 

(9) 

( 10 ) 


The pressure difference between both membrane sides 
causes a solvent velocity U inside the pore, which can be 
defined by the Hagen-Poiseuille equation (2). This 
assumption was validated elsewhere [9], [13]. 

U = r p 2 AP e / (Sr/Ax) (2) 

where AP e is the effective pressure (Nm' 2 ) and r] is the solvent 
viscosity (Nsm' 2 ). 

The ENP equation is different from the Nerst-Planck 
equation because it considers hindrance factors that are 
important to correct the convective and diffusive transport for 
a solute confined in a pore [13]. The values of hindrance 
factors depend on the ratio of ionic Stokes radius to 
membrane pore radius, named 2, and defined by (3), where r, 
is the Stokes radius. Expressions for the calculation of 
hindrance factors were proposed by Bowen and Mohammad 
and are represented by (4) and (5) for 0 < Xi < 0.8 [14]. 


Differentiation of (9) with respect to v and multiplication of 
(7) by n and summation over all ions give (11), which 
describes the potential gradient [9]: 


dy/ 

dx 


Tl D, p 


K -Ei-V — 

‘• c RT ' r 2 


c, -C, 


p J 


F (RT)± Zi 2 


(ID 


The Donnan equilibrium theory relates the electrochemical 
potential in the bulk feed solution to that within the pores [9] 
generating a difference in ions distribution at both membrane 
sides. The simultaneous contribution of Donnan and steric 
effects is usually represented by Donnan-Steric partitioning 
equation at feed and permeate sides, where Ay/ D (0) and 
A\ff D (Ax) are the Donnan potential (V) at feed and permeate 
interfaces, respectively [9]: 
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r, c u v 




0 

Yi C, 


= O f exp 


HI j 


( 12 ) 

(13) 


where y, and y t ° are the activity coefficient within the pore and 
in the bulk and c t (0) and c t (Ax) are the concentration of ion i at 
pore entrance and outlet, respectively, in mol m' 3 . 

Steric partition coefficient cj), is defined as the relation 
between solute concentration at bulk solution and at 
membrane pore and can be related to the parameter defined 
inEq. (14): 


<D,=(l-4) 2 


(14) 


The dielectric exclusion is mainly attributed to alterations 
in the solvent dielectric constant due to the confinement of 
water molecules within the pore. The effect known as Born 
effect corresponds to the variation of solvation energy when 
an ion is transferred from the bulk solution to the membrane 
pores. When the dielectric constant of the solution confined 
inside the pores is lower than that of the bulk solution, the 
excess solvation energy is positive and the ions are rejected by 
the membrane pores [9], [17]. 

This phenomenon has been represented in different ways. 
Bowen & Welfoot [9] have considered that the orientation of 
the water molecules at pore walls would lead to a reduction in 
dielectric constant, creating an energy barrier to solvation of 
ions into the pores which would increase salt rejection. 

The DSPM-DE model [10] considers a dielectric exclusion 
term which is attributed to the interaction of the ion with the 
electrical charges, induced by the ion, at interface between 
materials of different dielectric constants (the membrane 
matrix and the solvent). These induced charges are called 
“image charges” [18]. However, this consideration was not 
included in this work since the small radius of NF pores 
makes the pore solvent dielectric approach that of membrane, 
reducing the effect of image forces while increasing the 
solvation energy barrier [9]. 

The consideration of dielectric exclusion due to a single 
layer of ordered water molecules was validated by Oatley et 
al. [19] for Desal-5-DK NF membrane and was included in 
the model as a solvation energy term which multiplies the 
right-hand side of (12) and (13). This gives (15) and (16) [9]: 
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(16) 


The solvation energy barrier A W t (J) is calculated from the 
Born model [9]: 


size is estimated as follows [9]: 



This equation was developed with geometrical arguments 
and taking bulk dielectric constant to be s b = 80 (water 
dielectric constant at 25°C). The dielectric constant of the 
layer of oriented water molecules s shall be determined 
experimentally. 

B. Membrane Charge Density 

When a membrane is in contact with an aqueous 
electrolytic solution, it acquires an electric charge by many 
possible mechanisms such as dissociation of functional 
groups and adsorption of ions from solution. Thus membrane 
charges are influenced by the type and the concentration of 
ionic species in an electrolytic solution. Those surface 
charges have an influence on the distribution of ions in the 
solution due to the requirement of the electroneutrality of the 
system so that there is an excess of counter-ions in the 
adjacent solution. This leads to the formation of an electric 
double layer [2], [20]. 

A potential difference is created between the surface and 
the solution due to the membrane charge. The electrical 
potential decreases towards the solution until electroneutrality 
is reached. The electric double layer thickness represented by 
Debye length X D (m) is defined as follows: 

A d =yRT £l) 6- h )/{2F 2 r) (19) 

where s 0 is the permittivity of free space (8.85419 x 10‘ 12 C 2 
J’ 1 m’ 1 ) and I is the ionic strength (mol m’ 3 ) defined by: 


/= 0 . 5 $>, 2 C, 


( 20 ) 


Membrane charge of polymeric membranes is generally 
negative at high pH values (above isoelectric points), is 
neutral at around pH 3 - 4 and switches to positive at low pH 
values (below isoelectric points). Furthermore the isoelectric 
points depend on the electrolyte concentration and the anions 
adsorption is prevalent on the membrane, since they show 
lower hydration radii than cations [2]. 

C. Osmotic Pressure Difference 

The effective pressure difference AP e is defined as the 
difference between the applied pressure AP and the osmotic 
pressure difference An: 

AP e =AP-A7r (21) 

The osmotic pressure difference An is calculated 
considering the feed and permeates concentrations and can be 
estimated with the van’t Hoff relation at low concentrations 

[15]: 


AW. = 


2 2 
Zi e 


%7T£ 0 r i 


K £ p 


'bj 


(17) 


The variation of the pore dielectric constant s p with pore 


Ax = RTY J {c i , f -C i J ( 22 ) 

Other models have been proven to successfully predict the 
osmotic pressure, as the Pitzer equation, which was found to 
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predict the osmotic pressure within experimental error from 
dilute solutions up to an ionic strength pf 6 M. However, the 
Pitzer equation requires many salt parameters that many times 
are unavailable [5]. The van’t Hoff equation is simpler and 
meets the requirement of having a predictive model with less 
dependence of experimental data. 

The osmotic pressure difference is considered in the 
Hagen-Poiseuille equation to correlate the permeate flux with 
the applied pressure (23) [9]. This equation was used to 
evaluate the effectiveness of the van’t Hoff relation to predict 
the experimental permeate flux. 


r p 2 

8/7 v Ax J 


(23) 


III. Methodology 

A. Membranes and Feed Solution 

The membranes considered in this work were NF90 (Dow 
FilmTech) and NP030 (Microdyn-Nadir). The main 
characteristics of these membranes are shown in Table I. 


Table I. Summary of NF90 and NP030 membranes 
characteristics 



NF90 

NP030 

Supplier 

Dow FilmTec 

Microdyn-N adir 

Material of skin 
layer 3 

Polyamide 

Polyethersulfone 

Max. temperature, 

°c a 

45 

95 

Pore radius, nm 

0.55 b 

0.93 b 

pH range 3 

3-9 

0-14 

a According to membranes supplier; b [5]; c [21] 


In order to investigate the model adequacy for the NF90 
polyamide membrane and for the NP030 poly ether sulfone 
membrane and to assess differences in exclusion mechanisms 
between them, theoretical ion rejections of diluted single salts 
solutions typically present in seawater (NaCl, MgS0 4 , CaS0 4 
and Na 2 S0 4 ) permeating through the NF90 and NP030 
nanofiltration membranes were compared to the experimental 
results of Nicolini et al. [22]. The same experimental 
conditions were used in the simulations. The feed 
concentration, pH and zeta potential measurements for each 
salt solution are showed in Table II. 


Table II. Permeation conditions for NF90 and NP030 
membranes 


Salts 

c f 

(mol L' 1 ) 

pH 

NF90 

zeta 

potential 

(mV) 

NP030 

zeta 

potential 

(mV) 

NaCl 

0.025 

5.83 

-33.06 

-47.73 

Na 2 S0 4 

0.025 

6.22 

-58.37 

-23.06 

MgS0 4 

0.025 

7.52 

-27.25 

-22.45 

CaS0 4 

0.025 

5.83 

-22.02 

-25.89 


The experimental performance of the nanofiltration 
membranes was evaluated using a membrane filtration set up 
described in [22]. The system was equipped with a feed tank 
(10 L) where the feed solution was kept at constant 
temperature (23°C) by using a thermostatic bath. Feed 


solution was circulated through the membrane cell in preset 
flow rates and operation pressures, adjusted by flowmeter and 
a needle valve, respectively. The tests were performed in 
recirculation mode to keep feed concentration approximately 
constant. 

A flatsheet membrane with an effective filtration surface of 
0.0028 m 2 was used. Membranes were cut into circular pieces 
and then soaked overnight in ultra-pure deionized water. 
Afterwards, they were compacted at maximum operation 
pressure with ultrapure water or electrolyte solution until 
constant permeate flow. Membranes permeabilities were 
determined with ultra-pure water, at a cross-flow rate of 40 
L.h" 1 , varying applied pressure from 5 to 20 bar. 

The effective membrane thickness, Ax, was directly 
obtained from the permeability of pure water (P m ) based on 
(24): 

Ax = —(24) 

8 r?P m 

Zeta potential values showed in Table II were determined 
by Nicolini et al. [22] by tangential streaming potential 
measurements with an electrokinetic analyzer (SurPASS, 
Anton Paar) with a clamping cell. For each measurement, 
membrane samples of 55 mm x 25 mm were mounted 
opposite each other and separated with a spacer. The flow of 
an electrolyte solution through the channel under pressure 
generated a streaming potential, which was used to calculate 
the zeta potential with the Fairbrother-Mastin equation, (25) 
[23]: 

Aff> _ - £p£ r f W| (25) 

a p ScAl r J 

where (is the zeta potential (V), A (pis the measured streaming 
potential in the flow cell (V), s r is the relative dielectric 
constant of the electrolyte solution, X 0 is the bulk conductivity 
of the circulating electrolyte (mS.m-1), r] 0 is the bulk solvent 
viscosity (N.s m' 2 ), R h and R are the measured electrical 
resistances (mV A' 1 ) across the flow channel filled with the 
saline reference solution and with the electrolyte solution, 
respectively. 


B. Concentration Polarization 

Rejections observed during the experiments are defined by 
observed rejection (R 0 bsY 

R h = \-SjL (26) 

obs ^ 

C / 

However, in the presence of concentration polarization, the 
actual concentration at the membrane entrance is the wall 
concentration (C w ), which is higher than the feed 
concentration (Cf). As a result, the real rejection of a solute 
(Rreai) is higher than the observed rejection and is defined as 
follows: 

R ,=1-Tl (27) 

real ^ 

^ w 

The wall concentration can be correlated to the feed 
concentration by (28): 


= exp 


r . a 

I 

\ k cP 


(28) 


r -C 
C -C 

where j v is the volumetric flux through the membrane and k c is 
the mass transfer coefficient in the polarized layer. 

Many works that employed cross-flow modules in the 
permeation experiments in the laboratory scales reported low 


9 


www.erpublication.org 















International Journal of Engineering and Technical Research (IJETR) 
ISSN: 2321-0869 (O) 2454-4698 (P), Volume-7, Issue-8, August 2017 


or negligible concentration polarization [2], [14]. In the 
present work, the mass transfer coefficient for the laminar 
cross-flow cell was calculated with (29) [24]. The real 
rejections will be used for comparison with calculated 
rejections. 

(29) 


Sh = = 1.85Re 0 33 Sc 033 


D, 


(<dJLT 


where d h is the hydraulic diameter (m), L is the channel length 
(m), Re is the Reynolds number, Sc is the Schmidt number 
and Sh is the Sherwood number. 

C. Model Description 

Fig.l illustrates the algorithm developed in this work for 
the simulation of the ion rejections using DSPM-based 
models. This flowchart represents the process pathway to 
predict the rejections and fluxes of salts solutions. The routine 
was implemented in the software Scilab and simulations 
results were compared to real experimental data in order to 
evaluate predictability of this model. 

The following assumptions are usually considered in 
unidimensional models [13] and were also considered in this 
work: (1) concentration gradient is considered only along the 
membrane thickness Ax; (2) the electroneutrality condition 
must be fulfilled in the feed solution, membrane and 
permeate; (3) ion partitioning between membrane and 
solution is determined by Donnan equilibrium, steric effects 
and dielectric effects; (4) charge density is uniform along the 
membrane and (5) solute and solvent transport take place in 
cylindrical pores of known effective radius r p . 

The first step for solving the models equations involves an 
initial knowledge of membrane, solvent and ions parameters 
as well as operating conditions like temperature and pressure. 
Resolution involves an iterative step as ions concentration 
profile inside membrane depends on permeate concentration 
due to boundary conditions, the value of which is not known. 

Water properties depend on temperature (7) and salinity 
(s). Equation (30) indicates the correlation [25] used in this 
work to estimate water viscosity in the bulk solution. 


1965 


?7 = 1.234x 10“ s expl 0.00212s + ^ 


(30) 


The properties of each of the ions involved in this study are 
summarized in Table III [5], [9]. These values were 
incorporated in the DSPM-based models. 

Once water, ions and membrane properties information are 
known, transport model parameters can be determined using 
equations shown previously in this work. 

Ions concentrations inside membrane pore entrance, c ( (0), 
can be calculated from feed concentrations, Q/, with 
partitioning equation due to Donnan equilibrium and steric 
effects. Ions concentrations are related to each other by 
electroneutrality. 

An initial guess of permeate concentration, C i>p , is 
necessary to solve the concentration profile inside the pore, 
represented by ENP equation. Thus ions concentrations inside 
membrane pore exit, cl Ax), can be determined. As cl Ax) is 
also in equilibrium with Q p , permeate concentration can be 
estimated from partitioning. This iteration procedure is 
repeated until convergence constraint is reached. At this 
point, it is possible to calculate the rejections and compare 
with experimental data. 



Fig. 1. Algorithm developed for the solution of models 

Table III. Stokes radius, diffusivities and partial molar 
volumes of ions [5, 9] 


Ions 

n 

(nm) 

Djoo 

(HT 9 m 2 s' 1 ) 

Vi 

(cm 3 mol' 1 ) 

Na + 

0.184 

1.333 

-1.20 

Ca 2+ 

0.310 

0.791 

-18.04 

Mg 2+ 

0.348 

0.720 

-21.57 

Cl 

0.121 

2.031 

17.82 

so 4 2 ' 

0.231 

1.062 

14.18 


The effective charge density model used in this work was 
determined by taking into account the solution ionic force, pH 
and zeta potential data. This model was chosen because it 
considers the pH information that has an important role in 
membrane charge formation, which is not taken into account 
in the rejection fitting models at different concentrations. This 
chosen model also depends less on experimental data and 
adjusting parameters than those which consider functional 
groups dissociation and ion adsorption terms. 

In this model it is assumed that the effective surface charge 
density of the membrane is similar to the surface charge 
density at the shear plane (g s ), which is calculated using the 
simplified Gouy-Chapmann equation [20]: 

£ o £ bC_ ( 31 ) 


cr. = 


Ap 


The membrane effective surface charge density can then be 
estimated using the zeta potential values. Assuming that the 
membrane surface charge is uniformly distributed in the 
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cylindrical pore, it is converted to concentration units by [12]: 
V, = {lcr s )/{r p F) (32) 

The zeta potential data obtained by measurements were 
used as input parameters to calculate membrane effective 
charge density. The simulations were carried on in the same 
conditions of experiments and the results were compared to 
investigate model accuracy. 

IV. Results and discussion 

A. Determination of Real Rejections 

The real rejections of salts were calculated based on the 
experimental observed rejections considering of the 
concentration polarization phenomena with (28). The results 
are presented in Table IV for maximum and minimum 
experimental pressures. 

Table IV. Comparison of real and observed rejections for 
NF90 and NP030 membranes 


NF90 NP030 

Salts AP - 


(bar) 

Robs 

Rreal 

AR 

Robs 

Rreal 

AR 


% 

% 

% 

% 

% 

% 

NaCl 20 

0.940 

0.961 

2.1 

0.299 

0.326 

2.8 

5 

0.932 

0.937 

0.5 

0.244 

0.250 

0.6 

MgS0 4 20 

0.988 

0.993 

0.5 

0.637 

0.683 

4.7 

5 

0.989 

0.990 

0.1 

0.550 

0.560 

1.0 

CaS0 4 20 

0.990 

0.995 

0.5 

0.695 

0.733 

3.9 

5 

0.995 

0.996 

0.1 

0.607 

0.618 

1.1 

Na 2 S0 4 20 

0.993 

0.996 

0.2 

0.787 

0.812 

2.5 

5 

0.992 

0.993 

0.1 

0.806 

0.812 

0.6 


It can be noted form Table IV that the effect of 
concentration polarization decreases with decreasing 
membrane flux, as expected, for both membranes and all salt 
solutions. These results also show that the concentration 
polarization was low for the studied conditions. The 
maximum deviations from real to observed rejections were 
found to be 2.1% for NF90 membrane and 4.7% for NP030. 

B. Permeate Flux Behavior 

The pure water permeability of NF90 and NP030 
membranes was 4.69 Lh'W^ar' 1 and 1.53 Lh'W^ar' 1 
respectively. The effective membrane thickness was 
determined from (29) and found to be 0.49 pm and 9.9 pm for 
NF90 and NP030 membranes respectively. The difference 
between the thicknesses of the two membranes reflects the 
low permeability of NP030 membrane as compared to the 
NF90 membrane. 

The calculated permeate fluxes were determined with (28). 
Fig. 2 shows the experimental and calculated permeate fluxes 
vs. applied pressure for each of the salts solutions and a 
comparison to the pure water flux through both membranes. 
The flux reduction is due to the osmotic effect which can be 
observed even at low concentrations. For all salts studied, 
there is a nearly linear relationship between calculated 
permeate fluxes and applied pressure, which was in 
agreement with the experimental data. Table VError! 
Reference source not found, shows the comparison between 
calculated and experimental values on the basis of average 


deviation for all salts studied. These results show 
experimentally that the use of Hagen-Poiseuille equation 
defined in (2) is reasonable and confirm that the osmotic 
pressure can be well estimated using van’t Hoff equation and 
should not be neglected. 

C. Salt Rejections 

Initial simulations were performed without considering the 
dielectric exclusion in order to investigate the importance of 
this phenomenon for each membrane. In these cases, pore 
dielectric constants s p were considered equal to bulk dielectric 
constant (s b = 80) for all salt solutions. 

For the simulations using the model with the dielectric 
exclusion consideration, the magnitude of dielectric constant 
of ordered water layer, s*, was reassessed for each single salt 
solution using experimental data. The salts permeation was 
modeled with different dielectric constants of the oriented 
water layer lying between the two limiting values (s =6, 
when the water inside the pore is completely polarized, and s* 
= 80, when there is no polarization) and the optimum values 
which gave the best fit for each salt solution were found to be 
33+2 and 67 + 5 for NF90 and NP030 membranes, 
respectively. Pore dielectric constants estimations were then 
calculated with Eq. (23) and are present in Table V. 
Comparison between calculated and experimental 
permeate fluxes on the basis of average deviations (Dev.) 


Salts 

AP 

bar 

NF90 



NP030 


Tr exp 

(10-6 

m/s) 

Jv,calc 

(10-6 

m/s) 

Dev 

(%) 

Tr exp 

(10-6 

m/s) 

Jv,calc 

(10-6 

m/s) 

Dev. 

(%) 

NaCl 

20 

24.73 

24.64 

0.37 

7.07 

6.95 

1.64 


15 

18.85 

18.13 

3.79 

5.23 

5.18 

1.02 


10 

12.37 

11.61 

6.15 

3.53 

3.44 

2.51 


5 

4.71 

5.09 

8.01 

1.77 

1.66 

5.76 

MgS04 20 

21.20 

21.15 

0.23 

7.53 

7.32 

2.77 


15 

17.67 

15.47 

12.44 

5.78 

5.43 

6.09 


10 

10.60 

9.96 

6.01 

3.26 

3.53 

8.53 


5 

4.71 

4.31 

8.54 

1.50 

1.65 

10.08 

CaS04 

20 

24.31 

23.21 

4.54 

7.07 

6.81 

3.69 


15 

20.12 

17.00 

15.54 

5.30 

5.04 

4.86 


10 

12.60 

10.93 

13.28 

3.53 

3.28 

7.19 


5 

6.23 

4.69 

24.72 

1.77 

1.51 

14.28 

Na2SO 

20 

21.20 

19.85 

6.38 

7.07 

6.52 

7.77 

4 

15 

16.49 

14.39 

12.74 

5.30 

4.78 

9.90 


10 

10.60 

8.96 

15.51 

3.53 

2.99 

15.52 


5 

4.71 

3.46 

26.51 

1.77 

1.23 

30.36 


Table VI. 

Fig. 3 shows the simulated rejections versus applied 
pressure for NaCl, Na 2 S0 4 , CaS0 4 and MgS0 4 solutions 
without and with the consideration of dielectric exclusion (no 
DE and with DE) for both NF90 and NP030 membranes. 
Regarding experimental data, Fig. 3 shows that although the 
solvent flux increases with pressure, the rejection remains 
practically constant for each salt solution. This suggests that 
ions fluxes also increase with pressure. The order of rejection 
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was Na 2 S0 4 > CaS0 4 = MgS0 4 > NaCl for both membranes, 
which is a consequence of anionic electrostatic repulsion and 
the preferential attraction of divalent cations. Increased 
divalent cations concentration reduces the electrical exclusion 


of anions and, therefore, the saline rejection. NaCl salt 
presented the lowest rejection due to a less anionic repulsion 
for a monovalent ion and lower hydrated ionic radii. 



20 - - 


Na2S04 Experimental / NF90 
Na2S04 Experimental / N P030 
Mode I w ith DE / N F90 
Pure water / N F90 
Model with DE/NP030 
Pure water / N P030 



(b) 


Applied Pressure (bar) 



(c) 


Applied Pressure (bar) 


(d) 


Applied Pressure (bar) 


Fig. 2. Experimental and calculated permeate fluxes as a function of applied pressure for NF90 and NP030 membranes: 

(a) NaCl, (b) Na 2 S0 4 , (c) CaS0 4 and (d) MgS0 4 




(a) 


Applied Pressure (bar) 


(b) 


Applied Pressure (bar) 


(c) 



(d) 



Fig. 3. Rejection of NF90 and NP030 membranes vs. applied pressure for different single salts solutions at 25 mmol.L 1 . The 
solid line represents rejection calculated with DE consideration and the dashed line represents rejection with no DE 
consideration: (a) NaCl, (b) Na 2 S0 4 , (c) CaS0 4 and (d) MgS0 4 . 
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Table V. Comparison between calculated and experimental 
permeate fluxes on the basis of average deviations (Dev.) 


Salts 

AP 

bar 

NF90 



NP030 




7y, exp 

(ia 6 

m/s) 

Jv,calc 

(ia 6 

m/s) 

Dev 

(%) 

7y, exp 

(ia 6 

m/s) 

Jv, calc 

(ia 6 

m/s) 

Dev. 

(%) 

NaCl 

20 

24.73 

24.64 

0.37 

7.07 

6.95 

1.64 


15 

18.85 

18.13 

3.79 

5.23 

5.18 

1.02 


10 

12.37 

11.61 

6.15 

3.53 

3.44 

2.51 


5 

4.71 

5.09 

8.01 

1.77 

1.66 

5.76 

MgS0 4 20 

21.20 

21.15 

0.23 

7.53 

7.32 

2.77 


15 

17.67 

15.47 

12.44 

5.78 

5.43 

6.09 


10 

10.60 

9.96 

6.01 

3.26 

3.53 

8.53 


5 

4.71 

4.31 

8.54 

1.50 

1.65 

10.08 

CaS0 4 

20 

24.31 

23.21 

4.54 

7.07 

6.81 

3.69 


15 

20.12 

17.00 

15.54 

5.30 

5.04 

4.86 


10 

12.60 

10.93 

13.28 

3.53 

3.28 

7.19 


5 

6.23 

4.69 

24.72 

1.77 

1.51 

14.28 

Na 2 SQ 4 

20 

21.20 

19.85 

6.38 

7.07 

6.52 

7.77 


15 

16.49 

14.39 

12.74 

5.30 

4.78 

9.90 


10 

10.60 

8.96 

15.51 

3.53 

2.99 

15.52 


5 

4.71 

3.46 

26.51 

1.77 

1.23 

30.36 

Table VI. Dielectric constants from modeling of single salts 
solutions with NF90 and NP030 membranes 

Salt 


NF90 

£* 

£ p 


NP030 

£* 

£ p 


NaCl 


30 

42 


73 

76 


Na 2 S0 4 


36 

47 


63 

71 


MgS0 4 


33 

44 


67 

73 


CaS0 4 


32 

44 


63 

71 



For all salts, the model with dielectric exclusion (s < 80) 
predicted an increased rejection and led to a closer agreement 
between simulated and experimental data when compared to 
the DSPM model (e* = 80), indicating that the steric rejection 
and electrostatic partitioning alone are not capable of 
describing the rejection behavior. That emphasizes the 
importance of considering this effect in simulations. 

The calculated values for s* (Table VI) are all in the range 
30 to 36 for NF90 membrane and 63 to 73 for NP030 
membrane. The resulting dielectric constants of the oriented 
layer were found to be 33 ± 2 and 67 ± 5 for NF90 and NP030 
membranes, respectively. The result for NF90 membrane is in 
agreement with that report by [9] and [18] who obtained 
values of s*= 34.5 ± 2.5, s*= 35.5 ±1.5 and s*= 31 for the 
NF270, NF99HF and Desal-5-DK membranes respectively, 
which are also a polyamide membranes. The dielectric 
constants for NP030 membrane are interesting results, as they 
are quite different from those reported in the literature for 
polyamide membranes. This divergence can be attributed, in 
part, both to the differences in the membranes active layer 
polymer, as NP030 is a polyethersulfone membrane and to the 
differences in pore size, as NF polyamide membranes have 
very narrow pores, up to twice as less than those of NP030. 
This suggests that the interaction between the solvent and the 


membrane and the reorganization of solvent molecules in 
polyethersulfone membranes is quite different from the 
polyamide ones and the dielectric exclusion mechanism is less 
important since the magnitude of pore dielectric constant 
remains closer to that of the bulk solution. The reduction in 
the dielectric constant for both membranes, even if in 
different magnitudes, also supports the hypothesis of a single 
layer of ordered water molecules with modified properties 
and suggests that the inclusion of the solvation energy term is 
suitable for describing the dielectric exclusion for both NF90 
and NP030 membranes. 

The deviations of average rejections predicted by the 
model from those obtained experimentally are present in 
Table VII. 


Table VII. Comparison of simulated and experimental 
average rejections for NF9Q and NP03Q membranes 


Salt 

Average Deviation (%) 

NF90 

NP030 

NaCl 

1.51 

16.77 

Na 2 S0 4 

0.13 

11.38 

CaS0 4 

0.19 

6.91 

MgS0 4 

0.03 

6.77 


Good agreement between simulated and experimental 
results was observed with the model which considers the 


dielectric exclusion for both membranes. The accuracy 
obtained in the prediction of the experimental results using 
this consideration ranged from approximately 98.5% (for 
NaCl) to 100% (for MgS0 4 ) for NF90 and from 
approximately 83.2% (for NaCl) to 93.2% (for MgS0 4 ) for 
NP030 membranes. In the cases of CaS0 4 and MgS0 4 
solutions for both membranes, deviations from the 
experimental rejections were less than in the cases of Na 2 S0 4 
and NaCl solutions. It should also be noted that the model was 
able to provide more accurate results for the polyamide 
membranes (NF90) than for polyethersulfone membranes 
(NP030). This is an information that engineers interested in 
predicting industrial performance of these commercial 
membranes should be aware of. 

V. CONCLUSIONS 

In the present work ion rejection by nanofiltration of four 
salt solutions was investigated by using a DSPM-based model 
with and without considering the dielectric exclusion 
mechanism. Zeta potential data was used to determine the 
membrane charge density, osmotic effects were considered 
through the introduction of van’t Hoff equation in the model 
and pore dielectric constants were adjusted with experimental 
data. The aim of this study was to investigate the 
predictability of separation of the simple salts solutions of 
ions typically found in seawater by comparing the simulated 
rejection results with observed experimental data for two 
nanofiltration membranes, NF90 (polyamide membrane) and 
NP030 (polyethersulfone membrane). The concentration 
polarization was evaluated and the deviations between real 
and observed rejection were found to be low for both 
membranes. The results of predicted rejection and flux 
showed a good agreement when considering the dielectric 
effect by introducing the Born term in the model for both 
membranes, especially for NF90 membrane. However, when 
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this effect was not taken into account, the model 
underestimated the rejections for all salts solutions. This 
highlights the importance of considering the dielectric 
exclusion in the model and confirms that the determination of 
membrane charge density with zeta potential data can be 
considered. The dielectric constant of the oriented water layer 
was determined for each solution and found to be s =33+2 
and s = 67 + 5 for NF90 and NP030 membranes respectively, 
supporting that Born model is suitable to describe the 
dielectric exclusion for these membranes. The permeate flux 
behavior was also compared on the basis of average 
deviations which was considered low for each of the salts 
solutions at all pressures showing that the osmotic effect 
cannot be neglected in the model and that the van’t Hoff 
equation can provide acceptable results. 
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